MIT-CTP-3617 
hep-th/0509050 

Gravitational Correction to Running of Gauge Couplings 

Sean P. RobinsorQ and Frank Wilczeld 
Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics, 
Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

We calculate the contribution of graviton exchange to the running of gauge couplings at lowest non- 
trivial order in perturbation theory. Including this contribution in a theory that features coupling 
constant unification does not upset this unification, but rather shifts the unification scale. When 
extrapolated formally, the gravitational correction renders all gauge couplings asymptotically free. 
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The straightforward framework for quantum gravity — 
general relativity quantized for small fluctuations around 
flat space — is a famously nonrenormalizable quantum 
field theory 0, 0, 0, Q. Nevertheless, this framework 
is appropriate for describing interactions at energies and 
momenta below the Planck scale Mp = \/ he/ G-Ncwton ~ 
1.4 x 10 19 GeV/c 2 when treated as an effective low-energy 
theory. Indeed, if one makes subtractions to normalize 
physical couplings at an energy scale Eq well below Mp 
in such a way as to enforce the Einstein- Hilbcrt action of 
general relativity at the classical level with minimal cou- 
plings and a vanishing (or very small) cosmological term, 
then quantum corrections to this classical action at scale 
E will occur with coefficients containing positive powers 
of (E — Eo)/Mp, a small number. That procedure is the 
implicit foundation for practical use of classical general 
relativity as a model of nature despite the existence of 
quantum mechanics. It therefore underlies an enormous 
range of successful physical and astrophysical applica- 
tions. Only the classical theory really comes into play in 
those applications, because the quantum corrections are 
quantitatively small. Thus, the conceptual framework 
of effective field theory provides a sophisticated rational- 
ization for proceeding naively in applying the classical 
theory. 

Still, as Donoghue has emphasized 0, calculating cor- 
rections to the classical theory is a problem of method- 
ological interest. Moreover, quantitative considerations 
concerning interactions at ultra-high energy scales, per- 
haps approaching the Planck scale, are important in as- 
sessing the possibility of gauge theory coupling unifica- 
tion 0,0. Also, the size of gravitational corrections, in 
comparison to the leading classical term, give an objec- 
tive indication for the characteristic scale for the onset 
of quantum gravity phenomenology. With these motiva- 
tions, we consider here the one-loop (that is, first non- 
trivial order in perturbation theory) gravitational correc- 
tion to running of gauge theory couplings. 

We will perform this calculation directly in the frame- 
work described above. Any would-bc fundamental theory 
of quantum gravity should reproduce the same result in 
the limit of the physical scenario considered here, which 
is bosonic gravity in a four dimensional Minkowski back- 




FIG. 1: A typical Feynman diagram for a gravitational pro- 
cess contributing to the renormalization of a gauge coupling 
at one-loop. Curly lines represent gluons. Double lines repre- 
sent gravitons. The three-gluon vertex ■ is proportional to g, 
while the giuon-graviton vertex • is proportional to E/Mp. 

ground, with general matter and gauge sectors, at en- 
ergies below the Plank scale. Related calculations have 
been done in string theory 00, but this brings in several 
additional structures simultaneously, and we have found 
the results difficult to compare. 

Form of the correction. — The character of the cor- 
rection can be determined on very general grounds. 
The one-loop Feynman diagrams of interest involve a 
gluon vertex dressed by graviton exchange (See Fig. 
Alternatively, one could calculate the running coupling of 
a gluon to a test "matter" field. Gauge invariance (i.e., 
universality of the gauge coupling) implies that the same 
result must be obtained. This consideration highlights a 
cancellation between vertex and wave function renormal- 
ization, guaranteed by Ward identities, as is familiar in 
QED. 

Since the gauge boson vertex has strength g and gravi- 
tons couple to energy-momentum with a dimensional 
coupling oc X/Mp, dimensional analysis implies that the 
running of couplings in four dimensions will be governed 
by a Callan-Symanzik /3 function of the form 
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where the first term includes the familiar non- 
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gravitational contribution, and the second term includes 
the gravitational contribution. Since gravitons carry no 
gauge charge, b has the same value it had in the absence 
of gravitation, as determined by the matter content and 
the gauge group. Detailed calculation is required to de- 
termine the numerical value of the coefficient do- 

Even before knowing the value of do, much phe- 
nomenology can be extracted from the form of Eq. 0J. 
The equation can be integrated to yield 
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This reduces, of course, to the familiar logarithmic run- 
ning of inverse couplings in the limit oq — > (or Mp — > 
00). Corrections to the familiar result will be very small 
for E,E < M P . 

The value of clq is manifestly independent of the gauge 
interaction involved. So, if we consider several gauge 
couplings g 4 , with different b l , the condition that they 
unify at a common value g(E\j) is 
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for all choices of i and j. In particular, the left-hand 
side of Eq. J3J must be independent of the choice of 
Since that combination of initial couplings and renormal- 
ization group coefficients is not affected by the gravita- 
tional correction, the unification constraints remain un- 
changed. However, the scale of unification and the value 
of the common coupling at unification do change, as we 
shall now discuss. 

Comparing the unification condition with and without 
the gravitational correction, and taking E 2 / Mp 2 — » 0, we 
see that the relationship between the uncorrected unifi- 
cation energy E* and the corrected unification energy E\j 
is 
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If E* -C Mp, the self-consistent correction is 
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In standard (quasiminimal) unification schemes we find 
that E* is smaller than Mp by roughly 3 orders of mag- 
nitude, so this approximation is appropriate, but the cor- 
rection itself is of no practical importance. On the other 
hand, it is widely viewed as disturbing to have the sepa- 
ration of scales E* <C Mp . Theories that address this is- 
sue will inevitably bring in the gravitational correction — 
which, if do happens to be negative, helps to close the 
gap. 



The value of the coupling at unification is modified 
according to 
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where g* is the running coupling as determined by Eq. 
(JTJ with a -> 0. For E* -C M P , 
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Method of calculation. — The algebra required to eval- 
uate do in any straightforward way is formidable. We 
document algebraic details elsewhere 0b here we just 
sketch our method and conventions. 

The dynamics for a non-Abelian gauge field coupled 
to gravity in 3 + 1 spacetime dimensions is given by the 
action 
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where g = detg a f,, g a b is the spacetime metric, n 2 — 
16ir/Mp 2 , R is the Ricci scalar, g is the gauge coupling, 
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is the field strength, A* is the gauge field, / abc are the 
structure constants of the non-Abelian gauge group G, 
and A a is the spacetime covariant derivative operator. 
Since J-* b is antisymmetric under a «-> 6, the Christoffel 
connections arising from the derivatives in Equation @ 
cancel against each other, so the covariant derivatives 
here can be replaced with ordinary derivatives. 

We apply the background field method because it is 
especially well suited to the specific problem at hand in 
ways that we will highlight throughout this section. In 
accordance with this method, we seek to evaluate the 
effective action for classical field configurations by inte- 
grating over quantum fluctuations h ab and A^: 
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Here, we have expanded g a fc as a quantum fluctuation 
h ab about a background g a b, 



Sab — gab + hab, 
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and likewise expanded A* as a fluctuation A* about a 
background 
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In principle, the classical fields g ao and could satisfy 
the classical equations of motion, the coupled Einstein- 
Yang-Mills equations. For our purposes, however — that 
is, calculating the renormalization of gauge couplings to 
one-loop order in perturbation theory — it suffices to set 
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g a b equal to the flat Minkowski metric while allowing a a 
to obey the flat-space Yang-Mills equations of motion. 
We expand the action (JSJ in terms of these backgrounds 
and fluctuations up to quadratic order in the quantum 
fields, since within the background field method, higher- 
order terms in the action will only contribute to higher- 
loop processes. 

If we couple any matter to this system, we do not 
expand the matter fields as fluctuations about a back- 
ground because we are only interested in the renormal- 
ization of the gauge coupling. We still keep terms up to 
only quadratic order in quantum fields, however. Since 
we are expanding about Minkowski spacetime, matter 
terms in the action will give exactly the same contribu- 
tion to the one-loop renormalized coupling that they did 
in the absence of gravitation. 

The action JSJ is invariant under diffcomorphisms 

S v h a b =d a rib + dbVa + d a vi c h c b + d b r) c h ca + n c d c h ab , 



6 r ,A*=A*d a r 1 c + r) c d c A», 
5 v al =a a d a if + ifd c a a ; 



and under gauge transformations of the group G 

S a A a a =D 
S a h ab =0. 
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(Here, D a = d a — ia a t a when acting on representation r, 
and indices are raised and lowered with the background 
Minkowski metric.) These gauge symmetries need to 
fixed before performing the functional integration l|10|) . 
We take the background-covariant gauge-fixing condi- 
tions 



G a (A) =D a A aa = 0, 
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C a (h, A) =C a (h) - ^-F aab A a = 0. 
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where 



C a {h) = d b h ab - \d a h (h = h a a ), 
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and F a b is the appropriate function of classical fields only. 
Equation (|16fl is similar to an gauge Here it 

is engineered to cancel unpleasant graviton-gluon cross- 
terms that would otherwise appear in the expansion. Us- 
ing the Faddeev-Popov method |l2l | in conjunction with 
Feynman-'t Hooft weighting factors, the gauge choices 
each add gauge-fixing terms to the action as well as 
corresponding ghost fields. The ghost fields will not 
be expanded about a background and always appear at 
quadratic order in the action. So, like matter fields, the 
ghost fields contribute to the renormalized gauge cou- 
plings exactly as they do in the absence of gravitation. 
In particular, this means that the diffeomorphism ghost 



does not contribute at all and can be ignored in this back- 
ground field calculation. 

Not all gluon-graviton cross terms in the action can 
be eliminated by the choice of gauge (|16[) because gluon- 
graviton mixing in a vector background is a real physi- 
cal effect. In order to evaluate the Gaussian integrals in 
Eq. (|10f) as functional determinants, we formally combine 
h a b and A a into a "superfield" such that the mixing terms 
appear in the off-diagonal entries of the functional matrix 
in question. If multiple gauge symmetries are present, 
each with its own gauge field and renormalizable cou- 
pling, the superfield must be expanded to include each 
gluon type, as well as the graviton. The functional ma- 
trix then contains cross terms that mix different gluon 
types, but these do not ultimately contribute to the cal- 
culation at one-loop order. So, to this order, each gauge 
coupling gets renormalized independently. 

Result. — At this point the Gaussian integrals over the 
quantum fields in Eq. (|10J) are formally defined, but the 
resulting functional determinants contain ultraviolet di- 
vergences. We subtract them at a reference energy Eq. 
We find the one-loop effective action at energy scale E is 
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where bo depends on the gauge and matter content inde- 
pendently of whether gravitation is included in the cal- 
culation. Taking E differentially close to Eq , we read off 
the one-loop (3 function 
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Using k 2 = 16n/Mp 2 , the unknown coefficient in Equa- 
tion Q is now determined to be o,q = — 3/ir w —0.95. 

Comments. — The magnitude \oq\ ~ 1 indicates that 
Mp = ^/ 'he 3 1 'GNewton does indeed give a fair estimate of 
the energy scale for onset of quantum gravity, with no 
large numerical factors, for the problem considered here. 

At energy scales a few orders of magnitude below Mp 
the discussion that led to Eqs. JSJ) and @ is valid, so 
the negative sign of ao slightly increases E\j and slightly 
weakens the value of the unified coupling. This helps to 
close the gap between the unification scale and the Plank 
scale. 

Gravitational corrections will cause gauge couplings to 
run even in theories that in themselves are exactly con- 
formal invariant, that is when bo = 0. Two notable exam- 
ples in four dimensions are pure £7(1) electromagnetism 
and N = 4 Super- Yang-Mills [H Oj, Hi For these the- 
ories, the exponential integral in Equation J2J) has zero 
coefficient, so we are left with 
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FIG. 2: When gravity is ignored, the three gauge cou- 
plings of a model theory evolve as the inverse logarithm of 
E at one-loop order (dashed curves). Initial values at 100 
GeV were set so that the curves exactly intersect at approxi- 
mately 10 16 GeV. When gravity is included at one-loop (solid 
curves), the couplings remain unified near 10 16 GeV, but 
evolve rapidly towards weaker coupling at high E. 

That is, the coupling runs down from its infrared value 
as a Gaussian with a width of order Afp. For the pure 
U{\) case, this Gaussian running represents the renor- 
malized coupling strength of photons to non-dynamical 
or heavy sources, and would be the dominant — but still 
negligible — source of running in QED far below electron- 
positron threshold. For the theoretical "application" of 
Eq. J3U|) to N = 4 Super- Yang-Mills, and in the context 
of unification, it would be logical to include the contri- 
bution of gravitino-gluino loops to oq, but we have not 
calculated that here. 

The negative sign of ao also signifies that the gravi- 
tational correction works in the direction of asymptotic 
freedom: it causes the couplings to decrease at large en- 
ergy. Of course, its effect only becomes quantitatively 
important when the energy approaches the Planck scale, 
and soon after that one loses the right to neglect higher- 
order graviton exchanges. Though neglect of additional 
corrections is not justified beyond E <§C Mp, it is enter- 
taining to consider some consequences of extrapolating 
Eq. (J2J as it stands to these energies, taking into account 
ao < 0. The integral on the right-hand side converges 
as E — ► oo, and so Eqn. (|2(J[I arises as an asymptotic 
relation. Thus, the effective coupling vanishes rapidly 



beyond the Planck scale, rendering the gauge sector ap- 
proximately free at these energies. In Fig. O we illus- 
trate some aspects of the preceding discussion pictori- 
ally for an example theory with three gauge couplings 
whose low-energy values are chosen such that the b l de- 
termined from the matter sector result in a unification 
at E\j = 10 16 GeV. Obviously such a theory mimics the 
minimally supersymmetric standard model. 

SPR thanks Ted Baltz for useful discussions. This 
work is supported in part by funds provided by the U.S. 
Department of Energy (D.O.E.) under cooperative re- 
search agreement DE-FC02-94ER40818. 



* Electronic address: spatrick@mit.edu 
' Electronic address: wilczek@mit.edu 
[1] G. 't Hooft and M. Veltman, Annales Poincare Phys. 

Theor. A 20, 69 (1974). 
[2] S. Deser and P. van Nieuwenhuizen, Phys. Rev. D 10, 
401 (1974). 

[3] S. Deser and P. van Nieuwenhuizen, Phys. Rev. D 10, 
411 (1974). 

[4] S. Deser, H. S. Tsao, and P. van Nieuwenhuizen, Phys. 

Rev. D 10, 3337 (1974). 
[5] J. F. Donoghue, Phys. Rev. D 50, 3874 (1994). 
[6] H. Georgi, H. Quinn, and S. Weinberg, Phys. Rev. Lett. 

33, 451 (1974). 
[7] S. Dimopoulos, S. Raby, and F. Wilczek, Phys. Rev. D 

24, 1681 (1981). 
[8] E. Kiritsis and C. Kounnas, Nucl. Phys. B442, 472 

(1995). 

[9] E. Kiritsis, Introduction to super-string theory, Leuven 
notes in mathematical and theoretical physics, B9 (Leu- 
ven University Press, Leuven, Belgium, 1998). 

[10] S. P. Robinson, Ph.D. Thesis, M.I.T., 2005 (available 
at http://web.mit.edu/spatrick/www/PhDThesis/ 1; 
S. P. Robinson and F. Wilczek (to be published). 

[11] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B44, 189 
(1972). 

[12] V. N. Popov and L. D. Faddeev, BITP Report No. ITF- 
67-036, 1967 (unpublished); also Fermilab Report No. 
FERMILAB-PUB-72-057-T. 

[13] S. Mandelstam, Nucl. Phys. B213, 149 (1983). 

[14] M. T. Grisaru, M. Rocek, and W. Siegel, Phys. Rev. Lett. 
45, 1063 (1980). 

[15] D. R. T. Jones, Phys. Lett. 72B, 199 (1977). 



